SF3A2 promotes progression and cisplatin resistance in triple-negative breast cancer via alternative splicing of MKRN1.

Triple-negative breast cancer (TNBC) is the deadliest subtype of breast cancer owing to the lack of effective therapeutic targets. Splicing factor 3a subunit 2 (SF3A2), a poorly defined splicing factor, was notably elevated in TNBC tissues and promoted TNBC progression, as confirmed by cell proliferation, colony formation, transwell migration, and invasion assays. Mechanistic investigations revealed that E3 ubiquitin-protein ligase UBR5 promoted the ubiquitination-dependent degradation of SF3A2, which in turn regulated UBR5, thus forming a feedback loop to balance these two oncoproteins. Moreover, SF3A2 accelerated TNBC progression by, at least in part, specifically regulating the alternative splicing of makorin ring finger protein 1 (MKRN1) and promoting the expression of the dominant and oncogenic isoform, MKRN1-T1. Furthermore, SF3A2 participated in the regulation of both extrinsic and intrinsic apoptosis, leading to cisplatin resistance in TNBC cells. Collectively, these findings reveal a previously unknown role of SF3A2 in TNBC progression and cisplatin resistance, highlighting SF3A2 as a potential therapeutic target for patients with TNBC.

C9orf142 transcriptionally activates MTBP to drive progression and resistance to CDK4/6 inhibitor in triple-negative breast cancer.

BACKGROUND: Triple-negative breast cancer (TNBC) presents the most challenging subtype of all breast cancers because of its aggressive clinical phenotypes and absence of viable therapy targets. In order to identify effective molecular targets for treating patients with TNBC, we conducted an integration analysis of our recently published TNBC dataset of quantitative proteomics and RNA-Sequencing, and found the abnormal upregulation of chromosome 9 open reading frame 142 (C9orf142) in TNBC. However, the functional roles of C9orf142 in TNBC are unclear. METHODS: In vitro and in vivo functional experiments were performed to assess potential roles of C9orf142 in TNBC. Immunoblotting, real-time quantitative polymerase chain reaction (RT-qPCR), and immunofluorescent staining were used to investigate the expression levels of C9orf142 and its downstream molecules. The molecular mechanisms underlying C9orf142-regulated mouse double minute 2 (MDM2)-binding protein (MTBP) were determined by chromatin immunoprecipitation (ChIP) and dual-luciferase reporter assays. RESULTS: In TNBC tissues and metastatic lymph nodes, we observed that C9orf142 exhibited abnormal up-regulation, and its elevated expression was indicative of unfavorable prognosis for TNBC patients. Both in vitro and in vivo functional experiments demonstrated that C9orf142 accelerated TNBC growth and metastasis. Further mechanism exploration revealed that C9orf142 transcriptionally activated MTBP, thereby regulating its downstream MDM2/p53/p21 signaling axis and the transition of cell cycle from G1 to S phase. Functional rescue experiment demonstrated that knockdown of MTBP attenuated C9orf142-mediated tumour growth and metastasis. Furthermore, depletion of C9orf142 remarkably increased the responsiveness of TNBC cells to CDK4/6 inhibitor abemaciclib. CONCLUSIONS: Together, these findings unveil a previously unrecognized effect of C9orf142 in TNBC progression and responsiveness to CDK4/6 inhibitor, and emphasize C9orf142 as a promising intervention target for TNBC treatment.

ETHE1 Accelerates Triple-Negative Breast Cancer Metastasis by Activating GCN2/eIF2α/ATF4 Signaling.

Triple-negative breast cancer (TNBC) is the most fatal subtype of breast cancer; however, effective treatment strategies for TNBC are lacking. Therefore, it is important to explore the mechanism of TNBC metastasis and identify its therapeutic targets. Dysregulation of ETHE1 leads to ethylmalonic encephalopathy in humans; however, the role of ETHE1 in TNBC remains elusive. Stable cell lines with ETHE1 overexpression or knockdown were constructed to explore the biological functions of ETHE1 during TNBC progression in vitro and in vivo. Mass spectrometry was used to analyze the molecular mechanism through which ETHE1 functions in TNBC progression. ETHE1 had no impact on TNBC cell proliferation and xenograft tumor growth but promoted TNBC cell migration and invasion in vitro and lung metastasis in vivo. The effect of ETHE1 on TNBC cell migratory potential was independent of its enzymatic activity. Mechanistic investigations revealed that ETHE1 interacted with eIF2α and enhanced its phosphorylation by promoting the interaction between eIF2α and GCN2. Phosphorylated eIF2α in turn upregulated the expression of ATF4, a transcriptional activator of genes involved in cell migration and tumor metastasis. Notably, inhibition of eIF2α phosphorylation through ISRIB or ATF4 knockdown partially abolished the tumor-promoting effect of ETHE1 overexpression. ETHE1 has a functional and mechanistic role in TNBC metastasis and offers a new therapeutic strategy for targeting ETHE1-propelled TNBC using ISRIB.

MYC-driven U2SURP regulates alternative splicing of SAT1 to promote triple-negative breast cancer progression.

Triple-negative breast cancer (TNBC), although highly lethal, lacks validated therapeutic targets. Here, we report that U2 snRNP-associated SURP motif-containing protein (U2SURP), a poorly defined member of the serine/arginine rich protein family, was significantly upregulated in TNBC tissues, and its high expression was associated with poor prognosis of TNBC patients. MYC, a frequently amplified oncogene in TNBC tissues, enhanced U2SURP translation through an eIF3D (eukaryotic translation initiation factor 3 subunit D)-dependent mechanism, resulting in the accumulation of U2SURP in TNBC tissues. Functional assays revealed that U2SURP played an important role in facilitating tumorigenesis and metastasis of TNBC cells both in vitro and in vivo. Intriguingly, U2SURP had no significant effects on proliferative, migratory, and invasive potential of normal mammary epithelial cells. Furthermore, we found that U2SURP promoted alternative splicing of spermidine/spermine N1-acetyltransferase 1 (SAT1) pre-mRNA by removal of intron 3, resulting in an increase in the stability of SAT1 mRNA and subsequent protein expression levels. Importantly, spliced SAT1 promoted the oncogenic properties of TNBC cells, and re-expression of SAT1 in U2SURP-depleted cells partially rescued the impaired malignant phenotypes of TNBC cells caused by U2SURP knockdown both in vitro and in mice. Collectively, these findings reveal previously unknown functional and mechanism roles of the MYC-U2SURP-SAT1 signaling axis in TNBC progression and highlight U2SURP as a potential therapy target for TNBC.

Destabilization of microrchidia family CW-type zinc finger 2 via the cyclin-dependent kinase 1-chaperone-mediated autophagy pathway promotes mitotic arrest and enhances cancer cellular sensitivity to microtubule-targeting agents.

BACKGROUND: Microtubule-targeing agents (MTAs), such as paclitaxel (PTX) and vincristine (VCR), kill cancer cells through activtion of the spindle assembly checkpoint (SAC) and induction of mitotic arrest, but the development of resistance poses significant clinical challenges. METHODS: Immunoblotting and RT-qPCR were used to investigate potential function and related mechanism of MORC2. Flow cytometry analyses were carried out to determine cell cycle distribution and apoptosis. The effect of MORC2 on cellular sensitivity to PTX and VCR was determined by immunoblotting, flow cytometry, and colony formation assays. Immunoprecipitation assays and immunofluorescent staining were utilized to investigate protein-protein interaction and protein co-localization. RESULTS: Here, we identified microrchidia family CW-type zinc finger 2 (MORC2), a poorly characterized oncoprotein, as a novel regulator of SAC activation, mitotic progression, and resistance of cancer cells to PTX and VCR. Mechanically, PTX and VCR activate cyclin-dependent kinase 1, which in turn induces MORC2 phosphorylation at threonine 717 (T717) and T733. Phosphorylated MORC2 enhances its interation with HSPA8 and LAMP2A, two essential components of the chaperone-mediated autophagy (CMA) mechinery, resulting in its autophagic degradation. Degradation of MORC2 during mitosis leads to SAC activation through stabilizing anaphase promoting complex/cyclosome activator protein Cdc20 and facilitating mitotic checkpoint complex assembly, thus contributing to mitotic arrest induced by PTX and VCR. Notably, knockdown of MORC2 promotes mitotic arrest induced by PTX and VCR and enhances the sensitivity of cancer cells to PTX and VCR. CONCLUSIONS: Collectively, these findings unveil a previously unrecognized function and regulatory mechanism of MORC2 in mitotic progression and resistance of cancer cells to MTAs. These results also provide a new clue for developing combined treatmentstrategy by targeting MORC2 in combination with MTAs against human cancer.

Dynamic SUMOylation of MORC2 orchestrates chromatin remodelling and DNA repair in response to DNA damage and drives chemoresistance in breast cancer.

Rationale: SUMOylation regulates a plethora of biological processes, and its inhibitors are currently under investigation in clinical trials as anticancer agents. Thus, identifying new targets with site-specific SUMOylation and defining their biological functions will not only provide new mechanistic insights into the SUMOylation signaling but also open an avenue for developing new strategy for cancer therapy. MORC family CW-type zinc finger 2 (MORC2) is a newly identified chromatin-remodeling enzyme with an emerging role in the DNA damage response (DDR), but its regulatory mechanism remains enigmatic. Methods: In vivo and in vitro SUMOylation assays were used to determine the SUMOylation levels of MORC2. Overexpression and knockdown of SUMO-associated enzymes were used to detect their effects on MORC2 SUMOylation. The effect of dynamic MORC2 SUMOylation on the sensitivity of breast cancer cells to chemotherapeutic drugs was examined through in vitro and in vivo functional assays. Immunoprecipitation, GST pull-down, MNase, and chromatin segregation assays were used to explore the underlying mechanisms. Results: Here, we report that MORC2 is modified by small ubiquitin-like modifier 1 (SUMO1) and SUMO2/3 at lysine 767 (K767) in a SUMO-interacting motif dependent manner. MORC2 SUMOylation is induced by SUMO E3 ligase tripartite motif containing 28 (TRIM28) and reversed by deSUMOylase sentrin-specific protease 1 (SENP1). Intriguingly, SUMOylation of MORC2 is decreased at the early stage of DNA damage induced by chemotherapeutic drugs that attenuate the interaction of MORC2 with TRIM28. MORC2 deSUMOylation induces transient chromatin relaxation to enable efficient DNA repair. At the relatively late stage of DNA damage, MORC2 SUMOylation is restored, and SUMOylated MORC2 interacts with protein kinase CSK21 (casein kinase II subunit alpha), which in turn phosphorylates DNA-PKcs (DNA-dependent protein kinase catalytic subunit), thus promoting DNA repair. Notably, expression of a SUMOylation-deficient mutant MORC2 or administration of SUMO inhibitor enhances the sensitivity of breast cancer cells to DNA-damaging chemotherapeutic drugs. Conclusions: Collectively, these findings uncover a novel regulatory mechanism of MORC2 by SUMOylation and reveal the intricate dynamics of MORC2 SUMOylation important for proper DDR. We also propose a promising strategy to sensitize MORC2-driven breast tumors to chemotherapeutic drugs by inhibition of the SUMO pathway.

Protein Phosphatase 1 Subunit PPP1R14B Stabilizes STMN1 to Promote Progression and Paclitaxel Resistance in Triple-Negative Breast Cancer.

Triple-negative breast cancer (TNBC) represents the most lethal subtype of breast cancer due to its aggressive clinical features and the lack of effective therapeutic targets. To identify novel approaches for targeting TNBC, we examined the role of protein phosphatases in TNBC progression and chemoresistance. Protein phosphatase 1 regulatory subunit 14B (PPP1R14B), a poorly defined member of the protein phosphatase 1 regulatory subunits, was aberrantly upregulated in TNBC tissues and predicted poor prognosis. PPP1R14B was degraded mainly through the ubiquitin-proteasome pathway. RPS27A recruited deubiquitinase USP9X to deubiquitinate and stabilize PPP1R14B, resulting in overexpression of PPP1R14B in TNBC tissues. Gain- and loss-of-function assays demonstrated that PPP1R14B promoted TNBC cell proliferation, colony formation, migration, invasion, and resistance to paclitaxel in vitro. PPP1R14B also induced xenograft tumor growth, lung metastasis, and paclitaxel resistance in vivo. Mechanistic investigations revealed that PPP1R14B maintained phosphorylation and stability of oncoprotein stathmin 1 (STMN1), a microtubule-destabilizing phosphoprotein critically involved in cancer progression and paclitaxel resistance, which was dependent on PP1 catalytic subunits α and γ. Importantly, the tumor-suppressive effects of PPP1R14B deficiency could be partially rescued by ectopic expression of wild-type but not phosphorylation-deficient STMN1. Moreover, PPP1R14B decreased STMN1-mediated α-tubulin acetylation, microtubule stability, and promoted cell-cycle progression, leading to resistance of TNBC cells to paclitaxel. Collectively, these findings uncover a functional and mechanistic role of PPP1R14B in TNBC progression and paclitaxel resistance, indicating PPP1R14B is a potential therapeutic target for TNBC. SIGNIFICANCE: PPP1R14B upregulation induced by RPS27A/USP9X in TNBC increases STMN1 activity, leading to cancer progression and paclitaxel resistance.

Inositol monophosphatase 1 (IMPA1) promotes triple-negative breast cancer progression through regulating mTOR pathway and EMT process.

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer, which is characterized by high heterogeneity and metabolic dysregulation. Inositol monophosphatase 1(IMPA1) is critical for the metabolism of inositol, which has profound effects on gene expression and other biological processes. Here, we report for the first time that IMPA1 was upregulated in TNBC cell lines and tissues, and enhanced cell colony formation and proliferation in vitro and tumorigenicity in vivo. Additionally, IMPA1 promoted cell motility in vitro and metastatic lung colonization in vivo. Mechanistic investigations by transcriptome sequencing revealed that 4782 genes were differentially expressed between cells with IMPA1 knockdown and control cells. Among the differentially expressed genes after IMPA1 knockdown, five significantly altered genes were verified via qRT-PCR assays. Morerover, we found that the expression profile of those five targets as a gene set was significantly associated with IMPA1 status in TNBC cells. As this gene set was associated with mTOR pathway and epithelial-mesenchymal transition (EMT) process, we further confirmed that IMPA1 induced mTOR activity and EMT process, which at least in part contributed to IMPA1-induced TNBC progression. Collectively, our findings reveal a previously unrecognized role of IMPA1 in TNBC progression and identify IMPA1 as a potential target for TNBC therapy.

TOLLIP-mediated autophagic degradation pathway links the VCP-TMEM63A-DERL1 signaling axis to triple-negative breast cancer progression.

Triple-negative breast cancer (TNBC) is the most challenging breast cancer subtype to treat due to the lack of effective targeted therapies. Transmembrane (TMEM) proteins represent attractive drug targets for cancer therapy, but biological functions of most members of the TMEM family remain unknown. Here, we report for the first time that TMEM63A (transmembrane protein 63A), a poorly characterized TMEM protein with unknown functions in human cancer, functions as a novel oncogene to promote TNBC cell proliferation, migration, and invasion in vitro and xenograft tumor growth and lung metastasis in vivo. Mechanistic investigations revealed that TMEM63A localizes in endoplasmic reticulum (ER) and lysosome membranes, and interacts with VCP (valosin-containing protein) and its cofactor DERL1 (derlin 1). Furthermore, TMEM63A undergoes autophagy receptor TOLLIP-mediated autophagic degradation and is stabilized by VCP through blocking its lysosomal degradation. Strikingly, TMEM63A in turn stabilizes oncoprotein DERL1 through preventing TOLLIP-mediated autophagic degradation. Notably, pharmacological inhibition of VCP by CB-5083 or knockdown of DERL1 partially abolishes the oncogenic effects of TMEM63A on TNBC progression both in vitro and in vivo. Collectively, these findings uncover a previously unknown functional and mechanistic role for TMEM63A in TNBC progression and provide a new clue for targeting TMEM63A-driven TNBC tumors by using a VCP inhibitor.Abbreviations: ATG16L1, autophagy related 16 like 1; ATG5, autophagy related 5; ATP5F1B/ATP5B, ATP synthase F1 subunit beta; Baf-A1, bafilomycin A1; CALCOCO2/NDP52, calcium binding and coiled-coil domain 2; CANX, calnexin; DERL1, derlin 1; EGFR, epidermal growth factor receptor; ER, endoplasmic reticulum; ERAD, endoplasmic reticulum-associated degradation; HSPA8, heat shock protein family A (Hsp70) member 8; IP, immunoprecipitation; LAMP2A, lysosomal associated membrane protein 2; NBR1, NBR1 autophagy cargo receptor; OPTN, optineurin; RT-qPCR, reverse transcription-quantitative PCR; SQSTM1/p62, sequestosome 1; TAX1BP1, Tax1 binding protein 1; TMEM63A, transmembrane protein 63A; TNBC, triple-negative breast cancer; TOLLIP, toll interacting protein; VCP, valosin containing protein.

Bulk RNA Sequencing With Integrated Single-Cell RNA Sequencing Identifies BCL2A1 as a Potential Diagnostic and Prognostic Biomarker for Sepsis.

Background: Sepsis is one of the leading causes of morbidity and mortality worldwide in the intensive care unit (ICU). The prognosis of the disease strongly depends on rapid diagnosis and appropriate treatment. Thus, some new and accurate sepsis-related biomarkers are pressing needed and their efficiency should be carefully demonstrated. Methods: Differential expression analysis and weighted gene co-expression network analysis (WGCNA) were applied to detect sepsis and monocyte/macrophage-related genes. Least absolute shrinkage and selection operator (LASSO) and random forest regression analyses were used in combination to screen out prognostic genes. Single-cell RNA sequence profiling was utilized to further verify the expression of these genes on a single cell level. Receiver operating characteristic (ROC) curve and decision curve analysis (DCA) were also applied to verify the diagnostic value of the target biomarkers. Results: The intersections of the genes detected by differential expression and WGCNA analyses identified 141 overlapping candidate genes that were closely related to sepsis and macrophages. The LASSO and random forest regression analyses further screened out 17 prognostic genes. Single-cell RNA sequencing analysis detected that FCGR1A and BCL2A1 might be potential biomarkers for sepsis diagnosis and the diagnostic efficacy of BCL2A1 was further validated by ROC curve and DCA. Conclusions: It was revealed that BCL2A1 had good diagnostic and prognostic value for sepsis, and that it can be applied as a potential and novel biomarker for the management of the disease.

HSP90 N-terminal inhibitors target oncoprotein MORC2 for autophagic degradation and suppress MORC2-driven breast cancer progression.

AIMS: MORC family CW-type zinc finger 2 (MORC2), a GHKL-type ATPase, is aberrantly upregulated in multiple types of human tumors with profound effects on cancer aggressiveness, therapeutic resistance, and clinical outcome, thus making it an attractive drug target for anticancer therapy. However, the antagonists of MORC2 have not yet been documented. METHODS AND RESULTS: We report that MORC2 is a relatively stable protein, and the N-terminal homodimerization but not ATP binding and hydrolysis is crucial for its stability through immunoblotting analysis and Quantitative real-time PCR. The N-terminal but not C-terminal inhibitors of heat shock protein 90 (HSP90) destabilize MORC2 in multiple cancer cell lines, and strikingly, this process is independent on HSP90. Mechanistical investigations revealed that HSP90 N-terminal inhibitors disrupt MORC2 homodimer formation without affecting its ATPase activities, and promote its lysosomal degradation through the chaperone-mediated autophagy pathway. Consequently, HSP90 inhibitor 17-AAG effectively blocks the growth and metastatic potential of MORC2-expressing breast cancer cells both in vitro and in vivo, and these noted effects are not due to HSP90 inhibition. CONCLUSION: We uncover a previously unknown role for HSP90 N-terminal inhibitors in promoting MORC2 degradation in a HSP90-indepentent manner and support the potential application of these inhibitors for treating MORC2-overexpressing tumors, even those with low or absent HSP90 expression. These results also provide new clue for further design of novel small-molecule inhibitors of MORC2 for anticancer therapeutic application.

Chromatin complexes subunit BAP18 promotes triple-negative breast cancer progression through transcriptional activation of oncogene S100A9.

Triple-negative breast cancer (TNBC) is a highly lethal disease due to aggressive clinical phenotype and the lack of validated therapeutic targets. Our recent quantitative proteomic analysis of 90 cases of TNBC tissues and 72 cases of matched adjacent normal tissues revealed that the expression levels of BPTF-associated protein of 18 KDa (BAP18), a component of the MLL1 and NURF chromatin complexes, were upregulated in TNBC tissues relative to normal tissues. However, the biological function and the underlying mechanism of BAP18 in TNBC progression remain unexplored. Here, we report that BAP18 promoted TNBC cell proliferation, migration, and invasion in vitro and xenograft tumor growth and lung colonization in vivo. Mechanistic investigations revealed that S100 calcium-binding protein A9 (S100A9), a member of the S100 protein family that is frequently upregulated in breast tumors and acts as an oncogenic driver in breast cancer progression, was a downstream target gene of BAP18. BAP18 was recruited to histone H3 trimethylation at lysine 4 (H3K4me3)-marked promoter of S100A9 and enhanced its promoter activities. Notably, knockdown of BAP18 by short hairpin RNA in TNBC cells suppressed xenograft tumor growth in mice, the noted effect was partially reverted by re-expression of S100A9 in BAP18-depleted cells. Taken together, these results suggest that BAP18 promotes TNBC progression through, at least in part, transcriptional activation of oncogene S100A9, and represents a potential therapeutic target for TNBC.

RNF144A exerts tumor suppressor function in breast cancer through targeting YY1 for proteasomal degradation to downregulate GMFG expression.

Ring finger protein 144A (RNF144A), a poorly characterized member of the RING-in-between-RING family of E3 ubiquitin ligases, is an emerging tumor suppressor, but its underlying mechanism remains largely elusive. To address this issue, we used Affymetrix GeneChip Human Transcriptome Array 2.0 to profile gene expression in MDA-MB-231 cells stably expressing empty vector pCDH and Flag-RNF144A, and found that 128 genes were differentially expressed between pCDH- and RNF144A-expressing cells with fold change over 1.5. We further demonstrated that RNF144A negatively regulated the protein and mRNA levels of glial maturation factor γ (GMFG). Mechanistical investigations revealed that transcription factor YY1 transcriptionally activated GMFG expression, and RNF144A interacted with YY1 and promoted its ubiquitination-dependent degradation, thus blocking YY1-induced GMFG expression. Functional rescue assays showed that ectopic expression of RNF144A suppressed the proliferative, migratory, and invasive potential of breast cancer cells, and the noted effects were partially restored by re-expression of GMFG in RNF144A-overexpressing breast cancer cells. Collectively, these findings reveal that RNF144A negatively regulates GMFG expression by targeting YY1 for proteasomal degradation, thus inhibiting the proliferation, migration, and invasion of breast cancer cells.

Compelling Evidence Linking CD40 Gene With Graves' Disease in the Chinese Han Population.

Mutations in CD40 have been widely reported to be risk factors for Graves' disease (GD). The gene, along with its cognate ligand CD40L, may regulate pro-inflammatory and immune responses. Rs1883832, located at the -1 position of the Kozak sequence, is the most well-studied single nucleotide polymorphism (SNP) of CD40, and has been confirmed to predispose those with the alteration to GD, regardless of ethnicity. Our genome-wide association study (GWAS) indicated that several SNPs, including rs1883832 located within the vicinity of CD40 were associated with GD in the Han Chinese population. Aiming at identifying the most consequential SNP and its underlying pathogenic mechanism, we performed a two-stage refined study on 8,171 patients with GD and 7,906 controls, and found rs1883832 was the most significantly GD-associated SNP in the CD40 gene region (PCombined = 9.17×10-11, OR = 1.18). Through searching the cis-expression quantitative trait locus database and using quantitative RT-PCR, we further discovered that the rs1883832 genotype can influence CD40 gene transcription. Furthermore, we demonstrated that rs1883832 is a susceptibility locus for pTRAb+ GD patients. In conclusion, the current study provides robust evidence that rs1883832 can regulate CD40 gene expression and affect serum TRAb levels, which ultimately contributes to the development of GD.

Assessment of Molecular Subtypes in Thyrotoxic Periodic Paralysis and Graves Disease Among Chinese Han Adults: A Population-Based Genome-Wide Association Study.

Importance: Thyrotoxic periodic paralysis (TPP) is a potentially lethal complication of hyperthyroidism. However, only 1 specific susceptibility locus for TPP has been identified. Additional genetic determinants should be detected so that a prediction model can be constructed. Objective: To investigate the genetic architecture of TPP and distinguish TPP from Graves disease cohorts. Design, Setting, and Participants: This population-based case-control study used a 2-stage genome-wide association study to investigate the risk loci of TPP and weighted genetic risk score to construct a TPP prediction model with data from a Chinese Han population recruited in hospitals in China from March 2003 to December 2015. The analysis was conducted from November 2014 to August 2016. Main Outcomes and Measures: Loci specifically associated with TPP risk and those shared with Graves disease and prediction model of joint effects of TPP-specific loci. Results: A total of 537 patients with TPP (mean [SD] age, 35 [11] years; 458 male) 1519 patients with Graves disease and no history of TPP (mean [SD] age, 38 [13] years; 366 male), and 3249 healthy participants (mean [SD] age, 46 [10] years; 1648 male) were recruited from the Han population by hospitals throughout China. Two new TPP-specific susceptibility loci were identified: DCHS2 on 4q31.3 (rs1352714: odds ratio [OR], 1.58; 95% CI, 1.35-1.85; P = 1.24 × 10-8) and C11orf67 on 11q14.1 (rs2186564: OR, 1.50; 95% CI, 1.29-1.74; P = 2.80 × 10-7). One previously reported specific locus was confirmed on 17q24.3 near KCNJ2 (rs312729: OR, 2.08; 95% CI, 1.83-2.38; P = 8.02 × 10-29). Meanwhile, 2 risk loci (MHC and Xq21.1) were shared by Graves disease and TPP. After 2 years of treatment, the ratio of persistent thyrotropin receptor antibody positivity was higher in patients with TPP than in patients with Graves disease and no history of TPP (OR, 3.82; 95% CI, 2.04-7.16; P = 7.05 × 10-6). The prediction model using a weighted genetic risk score and 11 candidate TPP-specific single-nucleotide polymorphisms had an area under the curve of 0.80. Conclusions and Relevance: These findings provide evidence that TPP is a novel molecular subtype of Graves disease. The newly identified loci, along with other previously reported loci, demonstrate the growing complexity of the heritable contribution to TPP pathogenesis. A complete genetic architecture will be helpful to understand the pathophysiology of TPP, and a useful prediction model could prevent the onset of TPP.

Identification of BACH2 as a susceptibility gene for Graves' disease in the Chinese Han population based on a three-stage genome-wide association study.

The BACH2 gene regulates B cell differentiation and function and has been reported to be a shared susceptibility gene for several autoimmune diseases. Our previous genome-wide association study (GWAS) indicated that several single nucleotide polymorphisms (SNPs) in the BACH2 gene are associated with Graves' disease (GD) in the Chinese Han population; however, the association did not achieve genome-wide significance levels. Recently, this association of BACH2 with GD was confirmed in Caucasians in the UK population, but fine mapping in this region has not yet been reported. Here, we provide a refined analysis of a 331-kb region in the BACH2 gene, which harbors 359 SNPs, using GWAS data from 1,442 GD patients and 1,468 controls. The SNPs rs2474619 and rs9344996 were implied as the independent variants associated with GD by forward and two-locus logistic regression analysis. We genotyped eight out of 10 tagSNPs with P < 1 × 10(-3) in 3,508 GD patients and 3,209 controls, the results also showed that rs2474619 was independently associated with GD in the combined population from GWAS and the second stage (P = 1.81 × 10(-5)). The rs2474619 and rs9344996 were further genotyped in the third stage cohorts, and rs2474619 showed evidence of association with GD at genome-wide significance levels in the combined population (P = 3.28 × 10(-8), odds ratio = 1.13). The association of rs9344996 with GD can be explained by its linkage to rs2474619 in the combined population. Our study clearly demonstrated that BACH2 is a susceptibility gene for GD in the Chinese Han population and further supported rs2474619, in intron 2 of BACH2, is the best association signal with GD. However, the mechanism by which BACH2 confers increased risk of GD requires further study.

Refined association of TSH receptor susceptibility locus to Graves' disease in the Chinese Han population.

BACKGROUND: Convincing evidence has demonstrated the association of TSH receptor (TSHR) with Graves' disease (GD) in the Chinese Han population. OBJECTIVE: The aim of this study was to identify the causal variants for GD in the region encompassing TSHR by a refining association study. DESIGN AND METHODS: GD patients (1536) and 1516 sex-matched controls were recruited in the first stage, and an additional 3832 GD patients and 3426 sex-matched controls were recruited in the replication stage. Genotyping was performed using Illumina Human660-Quad BeadChips or TaqMan single nucleotide polymorphism (SNP) Genotyping Assays and the Fluidigm EP1 platform. RESULTS: When the results of regression analysis for 74 genotyped SNPs and 922 imputed SNPs in the first-stage cohort were combined, rs179243 and rs3783949 were the probable susceptibility SNPs associated with GD in TSHR. Eleven SNPs, including rs179243 and rs3783949, were selected to further refine the association in the replication study. Finally, rs12101261 and rs179243 were confirmed as independent GD susceptibility variants in the replication and combined populations. Further, we also found that the rate of persistent TSHR autoantibody positivity (pTRAb+) was significantly higher in the GD patients with the susceptible genotypes rs12101261 or rs179243 than in the GD patients carrying the protective genotypes, after the GD patients had been treated for more than 1 year. CONCLUSIONS: These findings indicate that rs12101261 and rs179243 are the possible causal SNPs for GD susceptibility in the TSHR gene and could serve as genetic markers to predict the outcome of pTRAb+ in GD patients.

Robust evidence for five new Graves' disease risk loci from a staged genome-wide association analysis.

Graves' disease (GD), characterized by autoantibodies targeting antigens specifically expressed in thyroid tissues causing hyperthyroidism, is triggered by a combination of genetic and environmental factors. However, only a few loci for GD risk were confirmed in the various ethnic groups, and additional genetic determinants have to be detected. In this study, we carried out a three-stage study in 9529 patients with GD and 9984 controls to identify new risk loci for GD and found genome-wide significant associations in the overall populations for five novel susceptibility loci: the GPR174-ITM2A at Xq21.1, C1QTNF6-RAC2 at 22q12.3-13.1, SLAMF6 at 1q23.2, ABO at 9q34.2 and an intergenic region harboring two non-coding RNAs at 14q32.2 and one previous indefinite locus, TG at 8q24.22 (Pcombined < 5 × 10(-8)). The genotypes of corresponding variants at 14q32.2 and 8q24.22 were correlated with the expression levels of C14orf64 and a TG transcript skipping exon 46, respectively. This study increased the number of GD loci with compelling evidence and indicated that non-coding RNAs might be potentially involved in the pathogenesis of GD.

A refined study of FCRL genes from a genome-wide association study for Graves' disease.

To pinpoint the exact location of the etiological variant/s present at 1q21.1 harboring FCRL1-5 and CD5L genes, we carried out a refined association study in the entire FCRL region in 1,536 patients with Graves' disease (GD) and 1,516 sex-matched controls by imputation analysis, logistic regression, and cis-eQTL analysis. Among 516 SNPs with P<0.05 in the initial GWAS scan, the strongest signals associated with GD and correlated to FCRL3 expression were located at a cluster of SNPs including rs7528684 and rs3761959. And the allele-specific effects for rs3761959 and rs7528684 on FCRL3 expression level revealed that the risk alleles A of rs3761959 and C of rs7528684 were correlated with the elevated expression level of FCRL3 whether in PBMCs or its subsets, especially in CD19(+) B cells and CD8(+) T subsets. Next, the combined analysis with 5,300 GD cases and 4,916 control individuals confirmed FCRL3 was a susceptibility gene of GD in Chinese Han populations, and rs3761959 and rs7528684 met the genome-wide association significance level (P(combined) = 2.27×10(-12) and 7.11×10(-13), respectively). Moreover, the haplotypes with the risk allele A of rs3761959 and risk allele C of rs7528684 were associated with GD risk. Finally, our epigenetic analysis suggested the disease-associated C allele of rs7528684 increased affinity for NF-KB transcription factor. Above data indicated that FCRL3 gene and its proxy SNP rs7528684 may be involved in the pathogenesis of GD by excessive inhibiting B cell receptor signaling and the impairment of suppressing function of Tregs.

Dense mapping of IL2RA shows no association with Graves' disease in Chinese Han population.

OBJECTIVE: Associations between IL2RA and various autoimmune diseases have been reported in Caucasians. We investigated whether genetic polymorphisms at the IL2RA locus were associated with Graves' disease (GD) in the Chinese Han population. DESIGN: We performed a genome-wide association study (GWAS) in 1 536 GD patients and 1 516 controls. The 1000 Genomes Project data were adopted as references for imputation analysis. After forward and conditional logistic regressions, we found that rs11256313 was the major risk variant in the CD25/IL2RA region. Thus, we further genotyped rs11256313 in a replication cohort with 3 694 GD patients and 3 510 controls using ABI 7900HT TaqMan Real-Time PCR System. RESULTS: Nine single nucleotide polymorphisms (SNPs) in the IL2RA block were nominally associated with GD in our GWAS (0·01 < P < 0·05). After imputation analysis, 13 imputed SNPs in the IL2RA block were weakly associated with GD (P ≤ 0·05). Logistic regression analysis suggested that the imputed rs11256313 could represent the IL2RA block (P = 0·003). However, we failed to replicate the association of rs11256313 in a larger cohort (P = 0·145). A subphenotype analysis of rs11256313 on thyroid hormone receptor antibody (TRAb) and gender showed that there was no association in any of the subphenotype groups (P > 0·05). CONCLUSIONS: The results suggested that common genetic polymorphisms at IL2RA do not exert a significant genetic effect on the development of GD in the Chinese Han population. Previously reported associations between CD25/IL2RA and autoimmune diseases including GD in Caucasians again imply that heterogeneity exists in different ethnic populations.